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ABSTRACT 


A general analysis of coplanar lines ona dielectric substrate is 
presented. The field equations are written in terms of Hertzian poten- 
tials and boundary conditions are applied to generate a system of 
linear equations. These equations are then Fourier transformed and 
reduced to a form where a computationally efficient solution can be 
obtained by the method of moments. Any desired accuracy may be 
obtained by expanding the transform of the current if the problem is 
one involving strip conductors or the electric field if there are slots. 
A one term expansion produces very good results. Both wavelength 
and characteristic impedance of the transmission line structure are 
obtained. Theoretical and experimental results are presented for 


coplanar strips. 
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I, INTRODUCTION 


During the past several years a variety of coplanar microwave 
transmission lines has been developed covering a wide spectrum of 
applications. However, there still remained several possible struc- 
tures to analyze, such as, coupled-slot lines, coplanar strips, etc. 
Also, theoretical analytic methods more accurate were sought " 
reduce the errors obtained by the more classical methods, such as 
the electrostatic approximation, the method of partial images, etc. 

The present work is concerned with two parallel, coplanar strips 
over a dielectric substrate, the other side being bare, as illustrated 
in Fig. 1. The structure is uniform and infinite in both x and z direc- 
tions; it is also assumed that the substrate material is lossless and 
its relative permittivity is € a. 

For coplanar strips (CPS) to be practical as a transmission line, 
radiation must be minimized. This is accomplished through the use 
of a high permittivity substrate, which causes the coplanar strip 
wavelength to be small compared to the free-space wavelength, and 
thereby results in the fields being closely confined to the strips with 
negligible radiation loss. 

A spectral domain transform method was snecatied by Professors 
Itoh and Mittra [Refs. 1, 2] which yields an exact solution; this method 
together with the method of moments [Ref. 3], is used in this present 


work. 
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Fig. 1. Parallel coplanar strips configuration 
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eeisePERSION* CHARACTERISTICS 


A. FIELD AND BOUNDARY CONDITIONS 
Let the following electric and magnetic fields exist along the 


propagation or z-direction, as justified by the Hertzian vector poten- 


tial functions described in Appendix A, 


dae Oo. - T2 


G2 C1) 
aus 2 p 2 gr 


@ 


(AZ) 


where ¥§ is the propagation constant, and since it was stated that the 


substrate material was lossless, it follows that 


a . (A3) 


From the field expressions in equations (Al) and (A2), all other 
components of electric and magnetic fields can be derived from 


Maxwell's curl equations, as shown in Appendix B, as, 
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Applying boundary conditions at the interface between regions 2 
and 3, tangential electric and magnetic fields must be continuous; 


therefore, aty = 0, 


Gs, (%,0,2) = Ge, (1,0,2) (A8) 
Gz, (%0,2) = Ex, (%,0,2) (A9) 
Wire 
Hz, (%,°, ) Ha, (v, 0,2) (A10) 
Hy, (2,0,2)= Hy, (4,0,2) 
(All) 


Similarly, at the interface between regions 1 and Z, tangential 
electric fields must be continuous and tangential magnetic fields 


discontinuous by corresponding surface current densities; therefore, 


at y = ai, 
Be, (¥,4,2) = Ga, (44,2) (A12) 
en @ae2 0624), ¥)*) 
. (A13) 
¥2 : 

da, (x,4,2)- de, (x,4,2) = Jy (x) e"~ om strips 

o elsewhere 
(A14) 

Ja (x) a on Strips 
Hy, (v,4,2)- Nx, (x,d, 2) = _ 


Oo elseurhere 


Also, the electric fields will exist only in the dielectric part of 


the interface and can be expressed as, 
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on Shibs 


(v,2,2) = (A 16) 
&2, ‘e. Gye™ elsewhere 

; © on strips 
Sr,(%, get e,(rye”* elsewhere Can) 


Substituting the field expressions of equations (Al) through (A7) 
into the boundary conditions expressions of equations (A8) through 
(A17) one obtains, 


Lo? & o° (x,0) = kes 6% G70) 












































(A118) 
ap® | (4 op oot (A) 
Y 2 (4,0) ~— |wp, = Eee) — | = L 
3x Jer ii lila acaaaiilan a alaarm by = 0) a9) 
(4) hl] 
he. Dae c 
2] o) = L.2 3, O'; iy ,0) (A20) 
ae . (e) (h) (e) 
Y = (X%,0) + WE, SCO) "= = (x,0} + JWE, 5 Src) (A21) 
ke? G(x, A) = k.? 5" (2, 3) (A2Z2) 
iy &) (e) (4) 
ab! 
(x, a) - JOpo = (¥,d) = = = (2,d)- {wpe — (x,¢) (A23) 
’ | 
fe, (ya)! - es iy oa) = Jy () (A24) 
ch) ¢e) 
y =e oer (x, 4) + J WEo oo (x, a) 
(A25) 
(4) (¢) 
|¥ 3 oe CY, d\ + Joer i (x, a) | = ae (x) 
hi? 5" (x, Aa) = Ca (x) . (A26) 
A 
y 2 iad op (A27) 


~ JeOPe rm = Gy (x) - 


1] 








B, SeweeCTRAL DOMAIN TRANS FORM 


All potential functions must satisfy the following relation: 


Wd+ke P=0 


(B1) 


where, 


, Tae ven Ep . a BR. 


(B2) 
One introduces the Fourier transform to the «-domain, as 


suggested by Itoh and Mittra [Refs. 1, 2], via, 


+ a9 . 
J®x , 
D; («,4) .f D; (x4) 2 dx ; t21,2,3 ; (B3) 
Therefore, one can find the tf-domain Fourier transform of 


equation (Bl) as, 


{2 FA [ir,y5] =- i, ise) ~3. ($5) 


B Ou® 
(B4) 


> ie &H(«,y) 
he? Bie, 4) = -(-jxl Bix) - —— 


oF Bimjy) = (0?~ ke?) Bay) soe 


Furthermore, one can define propagation constants for each 


region, 


25 7 + 2 7 
% 25 hei K +f ki (B6) 


where, 


de on WY PoE. 
k. - ki Ve 


(B7) 
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limemasplanar dielectric substrate theory, the effective permittivity 
is related to € - and € > as, 
Eo < Ege SS2 (B7a) 
Therefore, analyzing the expressions for propagation constants 


in equation (B6), for region 1 or 3, 


YP Ys = os Bie ky = wh cok (Fy (BS) 
where Us is the dielectric's effective wavelength and is related to 
the free-space wavelength A as, 

7 <b <a | (BB8a) 


Replacing this last expression into equation (B8), one can find 


2 
or ¥ Fas, 


lower and upper bound values for rf 1 


2 
oct 4 (e-4) (EP By? a ac? (B8b) 
therefore, § 4 and § 3 are always real quantities, independent of 
the values of & 


Similarly, for region 2, 


—_ 2 ~~ 
Yi = wry ge - ae, ( 5) ~ ene) 
: f 


where, as before, the upper and lower bound values for © are, 


<a yP es (1-¢,)( 2) 


(BY) 


; (B9a) 
lt is ekear that f 2 will be imaginary for small values of © and 


real for large values of & . Specifically, § , will be imaginary 


Z 
for 


x C oy €--1 (B9b) 


3 





-——-—_-—_——_— = —_— ee me — wm 
et A lil - es 


- 7 - 


- <_ aa 





anda real for 


Cid 
0 > = VErnd 


(Bee) 


Solving the potential functions differential equation (B5) for each 


of the three regions one obtains, as derived in Appendix C, 


field potential functions: 


ey =H Yea 
Oo" (Gy) = A ed @ 


6 Ck, if) = RB“ "Cmx) sinh Ry aE eu (ta) cosh fy , Y, = We 
al (x Re : (e) 
. (% 4) = Bro (x) sin Ye + C “(k) cos re >» Yerdxz 


(e) 
Os; (x, 4) = Dane e 
- Magnetic field potential functions: 


6" (x, 4 | AN sh 


if 


bh 
on (%, 4 ) Bg sinh Yay + o eS cosh hy, ¥, Re 


g" (h) 


DO, (uy, eS mc sin hy + CC) cos ey Ye= Dan 


im (4,4) = De @ ie 


where, clearly: 
¥, = Vs 
be, a kes 


and, 
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- Electric 


(B10) 


(B11) 


(Blla) 


(2) 


(B13) 


(B14) 


(B14a) 


(@3a5)) 


(B16) 


(Bl6a) 








BN) = | BCs) (B17) 
(e) _-¢ (2) 
B, (x)= yf Bc} ae, 
where both pi?) and ple) are evaluated at ¥ =j : 


In all the following discussion, equations subscripted with lower 


case "'a'' will refer to § , being imaginary. 


Z 
The set of equations (B10) through (B15) can be interpreted 

physically as follows: 

= Electric and magnetic fields are more concentrated in the 
dielectric substrate, region 2, and decay exponentially toward 
zero outwards into the surrounding media, regions 1] and 3. 

- Inthe dielectric substrate itself, electric and magnetic fields 
and, hence, energy, are more concentrated in the lower values 
of & . Although there is no direct comparison of the domains 

oO and xX , a set of plots of the potential functions indicate 
that energy is also confined to the vicinity of the strips. 
An approximate graphical representation of the energy distribu- 

tion within the dielectric substrate for a fixed value of dimension y 

is shown in Fig. 2. 

It is convenient to obtain the Fourier transform of equations 


(A118) through (A26) as follows, 


PS 





F(a) 


2x _-— 2k 
-= €--| = x €-~| 


Fig. 2. Energy distribution within the dielectric substrate 
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(B18) 
h.? BE! (w,0) = kes BS (x, 0} 





























ad, me (B19) 
- jx¥ Bz («,0) ~ Jwpo : (&,0) = ~jxY H.°'(%,0) ~ ,XPeo 34 (%, 0) 
b 
Lo? Bo) = ke? BY Ho) (B20) 
(e) ) 4) eese 
~ xr Gs" (x, 0) jeer = (e, ©) = ~jx¥ By (%,0} 7 or («,o} (B21) 
(2) 
b? oe Cr, a] a a D, (cg, (B22) 
(4) (4) 
4 = (B 3 
be? Bi cud) = ke? BE cma = Inco 
ie (B24) 
, (4) od 
-je ¥D, (&, ch) + J WEo 34 (x, A) 
(e} (B25) 
-[-juv BE o«,a) + joe 222 Ka) | = 
he. Bi cu, dA) = B> (ex) : (B26) 
. oh (el 
. (e) ' <2 = 
~ jx ¥ D, (x, A} - {pe a4 (%, A} &,. (x) an 
where the derivative transform pat, 
ob 
Pongal Je (B27a) 


has been applied. 
Upon substitution of the field expressions of equations (B10) 


through (B15) into the above relations, these become, 


Ld 





(B28) 
(é) (Z) 
ee Ce = ee? De) 


(B29) 
YOu) + Do Yo Bony = KY Ded + wy, % D™ (4 
B30 
12 Ge cosas ke? p") x) (B30) 
KY O(n) - wer %e Bay = KY Dy ~ We Ys Du) eee) 
£2 NN cx) = he? | Be sinh Yoda + AM, cosh Y¥p a | (EB S2) 
“KY A ces ~ Wo Y; Al”) cee) = & Y | a sinh %d + een cosh ea | 
(B33) 
y (h) (k) . , 
+ Oo Te E (x) cosh Yed + Cm) sinh ‘2 | 
b.2 ea Ke. Exc sinh Yoda + a (x) cosh Ye] = J (%) (B34) 
- jX¥ at) (0' ~ [Week Ae ce.) + jx [By sinh ted Hem Sh Yed 
(B35) 
~ JE Y> | BP en cosh te + Cn) sinh I, 3] = Ja (ew) 
he? Nay = Es (00) (B36) 
| 
; (e) (4) | 
JV ANT (ml + [wpe N A a S68) (B37) 


Equations (B28) through (B31) can be expressed in terms of 


eo and C (h) , obtaining: 
(c<) 





(c< ) 
2 
Dm) = ( <= | 0" Gee) (B38) 
(h) kee é (4) 
D (x) = ( <2) C (%) (B39) 
ke kez om 
Ben) me ome {er | 4- (<2 2 a | @ ea) + WE, Ya (= = yng » (B40) 
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sean for[t-(He}]cMensonns (EBC | cose 


J WE, Yo 


imi { r|( eee 1 | CF + Ope » (= ai oy | (B41) 


Be) x aa {uy [(G2)e4 : 1] ea + wpe Ya (=) Cea } (B4la) 
J oe 


where, in equations (B40a) and (B4la), 


No = 2 
. (e) (h) 
Furthermore, one can relate the constants C and C 
(O¢ ) (cz ) 
(e) (h) : 
interms of A (Q&)and A ‘(@&) via equations (B32) and B33), 
obtaining, 
B42 
aoc = PAM ot ~ & AM cout (ee6) 
oe (e! (k) 
CM Negy = Py AT ad = Ry ANC 00) (B43) 


(e) (h) 


A relation between the constants A §(o0)and A’ ‘() in terms 
of the surface current densities JL ) and J,() can be established 


via equations (B34) and (B35), obtaining, 


AM) « + [a Jo Ge) ~ Qy Sc Goo! | 
RQ2Qz + QiQy = en Bae 
(h) 1 | 
AX Je | 3 Jo («) + @ “i x) | 
aia,+ae. |” aw ee) 


Finally, applying equations (B36) and (B37), the surface current 
densities can be related to the electric fields & (oa) and & 2% ie 
x 


obtaining, 


‘ 
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er" eet” Gee 


ee VS 





— —_ = 
- em << 
_. tees -.— 
<-> a > —- eae 


| i ~« —_= —_— a m_— -—=_— = -~ 
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a” =p 








(B46) 


M, (%, 6) Jae) + Me (%, 8) Jx(d= Ga (a) 


B47 
Ms (<,6) Ja (x) + My (x, 8) ne Kis EC) ( ) 


where all the intermediate steps are fully derived in Appendix D. 


ee PrPHyYSsICAL PARAMETERS 

Up to this point, the problem has been stated ina quite general 
form, i.e., there is no dependence on the actual physical configura- 
tion other than the boundary surfaces. 

In general, the surface current densities, Jyh ) and J_(e< ), can 


be expressed as a train of known basis functions, [Ref. 3], as, 


Jy (x) 2 Li Ax; (ee) (C1) 
Ja Cm) = im Ge 42; Col . (C2) 


Substituting the above expressions into equations (B46) and (B47), 


one obtains, 


M, (%, 2) 2 b; 42; (el + No (ox, B) 2 i bx; (ul = Bs (a (C3) 
Ma (~, 6) 2 bi fa, (we) + My (%, 8) 2, Ai Ax; (&) = ey (ox) (C4) 


Since the two conductors are sufficiently narrow, one may assume 
that the surface current density in the x-direction is zero, which will 


modify equations (C3) and (C4) as, 
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M, («, 8) y bi Ai (x) = B. («\ (C5) 


cz] 


Ma (t, (3) 2s Gi fai (od = Ex (ul (C6) 


Furthermore, although theoretically one should perform a multi- 
term approximation to the surface current density in the z-direction, 
in this case a one term approximation is sufficient since a good 


assumption of this surface current density is a square pulse, i.e., 


Wr (%, 8) 5, fe Ce = Ba (0<\ (C7) 


Ma (~, 8) G, Ae, Col = &, C6} 
(C8) 


or, substituting equations (Cl) and (C2) for the case where i = l, 


By (%) (C9) 


By, (%) 


MM, (be, @) Ja (tc) 


Ns (%, 2) Js Ct) 
(C10) 


It is clear that one only needs to work with either one of the 
above equations; consider equation (C9), for the present case. 

As mentioned before, a good approximation to the z-direction 
surface current density is a square pulse; therefore, it is necessary 
to obtain the « -domain Fourier transform of JAX ), shown in 


PON, See uc 


aT (7A) 


Ztw 
Ja) = eran - J ot ay 
s 
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-8 
=~ Ww 


P| 





Nin 


tio 


hex 
ee 








1 Jjk® 
cen.“ 


a . & > S A 
[ett ee JXZ tee 


(C11) 


eee 
== 


Nit 


+W) 


0 


0 ac ~€2 


C. 


" 


s. A 
“i | 2 os x 2-2 cag < (S40) | 


J 
eee sin = (s+w) sin SW 
~ o< Z: 2 


Therefore, equation (C9) becomes, 


M, (2) Io gin & (sew sin Sw = Gym), (C12) 


In general, for a strip configuration, the electric fields are not 
easily defined; therefore, it is convenient to make equation (C12) 
independent of S (cc). This could be done by taking an inner product 
with a function ariacgemal to E(x ). Applying this concept, equation 


(C12) becomes, 


<M, (%, 8) Ja (ad, WCK) > = < Ga), Wie) > 


(G3) 
Define a suitable inner product, suchas, 
+ 00 
CV(K),WCK) > = J veel WICK) dn (C14) 
= 6O sf 


Ze, 





Pipe 3. 


+W 


Nin 


S 
=a —\VW = 


rin 
LS es) 


J,(*) 


Z - directed surface current density 
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A suitable weighting function W(t) could be the complex conjugate 
oun “ec ), i.e., J (- ¢). 

Z Fz 

Also, by Parseval's theorem, the right-hand side of equation 


(C13) becomes zero because of the orthogonality of the integrand. 


Therefore, equation (C13) becomes, 


+00 
te poe ss Ue = 
16 J My, (%, 3) aa SINS (S4+w) Sin 5S W dk =O | (C15) 


It is clear that once equation (C14) is integrated, the dependence 
on the variable ec disappears, so one can state that the whole process 
is really a function of frequency and the structure's physical charac- 
teristics, namely, width of the strips W, separation between the 
strips S, thickness of the substrate D, and the dielectric's relative 
Permittivity ¢ bas 

Once the e¢-domain Fourier transform of the surface current 
density is obtained, expressions for the different constants, A through 
D, have to be redefined to account for the imaginary character intro- 
duced by abe ). It turns out that all the magnetic potential constants, 
Peet een BY) ge), ae (eat and Tue are imaginary, and 
all the electric field constants, real. A detailed discussion of this 
procedure is presented in Appendix E. 

Using the expressions derived in Appendix E, the equations for 


the various constants are modified as, 





ery | a] 


C16 
[aal2,- aia) oe 
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—— eee 
rr ne —_ 
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> = EE eee . oa 












> ~~) = Eh Rie 


eo 2 ef anes = & 


2 
: -—— = —ai™ AX, 
=—<-— | Ga 
. [ie ; 
eanteliaed > =p 





(h 
ae 7 Raa (en 


la.les-@, Jax 
a ey 4 P AW) + |B | A! a (C18) 
Mee = 5 [IB] AM ~ Q [AM ca] (C19) 
ee (R) = = = 3 1 {wets (zy Oe) real ai ~ t]| cus] } (C20) 





eh =] : {«e[(E)- 1] a (oe) ONE: Tae ae “ces } (C21) 


Ge Ve 
DM cx} = (=) OO) Koes) 
Dewy = j (=) [oe | (C23) 


D. NUMERICAL INTEGRATION 

Essentially, the main problem, as far as numerical integration 
is concerned, is the evaluation of equation (C15). Although theo- 
metrically the limits of integration are - © and + o, truncation at 
points corresponding to sufficiently small values of the integrand is 
allowed. 

The numerical integration method used is a modified Simpson's 
rule. Consider a general function, as shown in Fig. 4, where one 
divides the axis of integration « in N segments of width A ; the 
ordinate, corresponding to the half-width point of each segment, is 


taken as the average ordinate of each integration segment and a 
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Fig. 4. Simpson's rule of numerical integration 
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Fig. 5. Typical integrand plot 
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rectangle is formed with the segment's width. It follows that the 
total area under the given curve will be approximately equal to the 
sum of the areas comprised by the N rectangles. 

It is clear that for more accurate results the width of each seg- 
ment A, should be a function of the rate of variation of the function 
E(&), i.e., for rapid-varying functions, the width AV should be 
smaller than for slow-varying functions. 

Several graphs of the integrand of the final equation (C15) were 
obtained; it is clear that for a given substrate of thickness D and 
relative permittivity € ,? one can vary the frequency of operation, 
width of the strip conductors W, and separation between them S, 
obtaining different corresponding values of effective wavelengths on 
the dielectric X . Therefore, there is no single graph that could 
be used throughout the complete analysis. However, a graph corres- 
ponding toa set of parameters chosen arbitrarily will give a fairly 
good idea of how the integrand behaves and, consequently, will 
indicate the approximate limits of integration and an acceptable 
segment's width A. A typical graph is shown in Fig. 5. 

It was chosen, according to the graphical expression of the inte- 
grand, that the limits of integration will be -1700 and +1700, since 
contributions beyond this point were less than 1% of the maximum 
contribution. Likewise, a good choice for the segment's width 


is 0.5, which will give a total of 6800 steps of integration. 
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The rate of change of the value of the integral in equation (C15) 
with respect to the effective dielectric wavelength IN was also 
investigated to determine whether it was a monotonically increasing 
or decreasing function, an oscillating function, or a non-uniform 
function. For this purpose, the integral was evaluated at several 
increasing values of X and it was observed that it behaved as a 
monotonically decreasing function, as shown in Fig. 6. 


The resulting integral in equation (C15), 


+1700 
{ - 2X ay ei OC 
{ M, (#8) ue SIN = (S+w) os W du 
~1760 


will reach the desired value of zero at an approximate value of 

for a given set of parameters, i.e., frequency, width of conductors, 
separation between conductors, thickness of the substrate, and 
dielectric permittivity. 


Since fs is defined as 
Ze 


> 
where XN varies as 
» 


ee ON SD 
Ver 


it is clear that an iteration procedure should be used, varying 


from a starting value equalto 2A q’ the dielectric's wavelength, 


nd oe 


Ver 


in appropriate steps toward A , until the value of zero is reached. 
The accuracy of the effective dielectric's wavelength depends 


critically on the size of the above mentioned iteration steps. 
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Therefore, there is a trade-off in practice of accuracy vice computer 
time or, equivalently, accuracy vice dollars. For the present work 
accuracies of the order of 0.1% were obtained by using iteration 
steps equal to 1/1000 of the corresponding free-space wavelengths. 
The integration subroutine itself was tried for several curves 
of known equations of various rates of variations, i.e., sinusoids, 
paraboloids, cubics, etc. In all cases variations of, at most, 0.2% 
their theoretical value were obtained. Also the limits of integration 
were increased to 2500 obtaining variations of, at most, 0.05%. 
Finally, integration steps of 0.4 and 0.6 were also tried to check 


accuracy and, again, variations of only 0.1% were observed. 


E. COMPUTER PROGRAMMING AND RESULTS 
1. Computer Program Organization 

The dispersion characteristic program, written in FORT- 
RAN IV language, accepts the following data: 

- Dielectric's permittivity 

- Width of the conductors, in millimeters 

. Thickness of the substrate, in millimeters 

- Ratio of the thickness of the substrate to free-space 

wavelength 
- Ratio of the separation between conductors to width of 


conductors. 


31 





Likewise, the limits of integration, the step of integration 
and the iteration step X, need to be specified. 

The program starts by calculating all the required physical 
constants, namely, NM, & , €, and po . All the «- and 


p -independent parameters are then evaluated, i.e., k Ko, a) , 


1? 
etc. The iteration step is initiated by assigning the lower bound 

value of A Mee NZ defining in this way the value of fs . At 
this point the integration is performed by assigning the lower bound 
value of & , i.e., -1700 for the present case defining consequently 


the parameters k.,, Kop, Kegs 6 4 ¥ >, and v x Means are 


1? 
provided to test the real or imagimary character of & 2 and two 
corresponding paths then follow. The integration, or better, the 
summation stops when & reaches its upper bound value, i.e., +1700 
for the present case, and means are provided to test whether the 
absolute value of the resulting quantity is less than an arbitrary 
small epsilon value, i.e., on for the present case. If the case is 
that the result is Tpeseee than this epsilon value, \ is increased 
eae yn in one iteration step X and the summation is performed 
again. The whole process ends either whena result less than epsilon 
is found, in which case the corresponding » is the answer sought, 


or when A reaches its upper bound value, i.e., XK 


The program is outlined in detail in Appendix F. 
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Dispersion characteristic curves for €= 12 
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Fig. 8. Dispersion characteristic curves for Se Jae 
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2. Computer results 
Effective dielectric wavelength ratio relative to free-space 
wavelength information was obtained for three different relative 
permittivities, namely, for €. = l2y hepwamdecOm hice info ramatwon is 
presented in terms of ratios of the different parameters S, W, D, and 
in order to permit the use of any set of actual physical parameters. 
For each permittivity, curves were obtained for the following 
values of the ratio of separation between conductors to width of con- 
ductors, S/W=:, 4, 1, 3/2, 2, 3, and4, corresponding to’values 
of the ratio of substrate thickness to free-space wavelength D/), 
perm 0.01 to 0.06. 
Dispersion curves for € sgl 12 are presented in Fig. 7, for 


Cc. = 16in Fig, 8, and for - =a Omi ip sees 


EF. EXPERIMENTAL RESULTS 

1. Fabrication 

The coplanar strip (CPS) fabrication was made by photo- 

etching technique. Photoresist is first spread uniformly across the 
surface and allowed to dry. This film is then exposed to ultraviolet 
light through a photographic mask that permits the light to fall only 
on those areas in which the copper layer is preserved. The un- 
exposed film is next washed away from the esas leaving a solid 
film of photoresist in the areas that were exposed. During the 


etching process the copper layer is removed in the unprotected areas 
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but is unaffected under the photoresist. Finally, the photoresist is 
removed with acetone, leaving only the dielectric surface with the two 
copper strips, as shown in Fig. 10. 

It was decided to construct the coplanar strips on a piece of 
6''x 3" x 1/8" copper clad substrate of relative permittivity of 12. 
Also, the ratio of width of conductors to thickness of the substrate 
was chosen arbitrarily as 1] and, similarly, the ratio of separation of 
conductors to width of conductors was also chosen arbitrarily as 1. 

The coplanar strips were fed with a transition consisting of 
a 2% inch piece of 85 mil semi-rigid coaxial line, soldered as shown 
aaecae. 10. 

2. Experimental data 

Basically, the measuring devices used were a slotted line 
fed by a sweep frequency oscillator and monitored by a voltage stand- 
ing wave ratio meter. 

From transmission line theory, it is known that an open -~ 
circuited line will have voltage nulls every half-wavelength, which 
in this case will be the effective dielectric wavelength. 

Using this procedure, data was collected for a range of 
frequencies from 1 to 6 GHz and a comparison was made with the 
corresponding computer results, observing an error of 0.1% which 
was, as a matter of fact, the accuracy that was claimed for these 


results. 
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Fig. 10. Coplanar strips fabrication 
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Table I lists the experimental results and compares them 


with the computer results. 
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TAS LE I 


DISPERSION CHARACTERISTIC COMPARATIVE RESULTS 


Coplanar strip of characteristics :€, = 12 


d 


1] 


OaZe 


s=w sd 


Frequency X/> from x7 from 
Giz Wales Fe experiment 
0 0.447 0.446 
Zt) 0.435 0.435 
30 0425 0.424 
4.0 0.415 | 0.414 
5.0 0.406 0.406 
6.0 O26 Oe 
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Ill, CHARACTERISTIC IMPEDANCE 


iio avERAGE POWER AS A FUNCTION OF THE HERTZIAN 
VECTOR POTENTIAL FUNCTIONS 


oF 


As stated by Collin [Ref. 4], a general expression for the time- 


average power flow in terms of the electric and magnetic fields is 


Re ny) = $@[ J Gxat- 2, as] (G1) 
where, 
(G2) 


& x" -G, = Cy = By Wr 
Referring back to equations (Al) through (A7), it is found that the 
equations for the electric and magnetic fields can be expressed in terms 


of Hertzian potential functions as, 














B= bE He” (G3) 
4 LE pee a“ oe 
., | (G4) 
ad@! ; ap" ¥2 
Bi ° (Y Ss Og Je (G5) 
age apm Y2 
Bq = (y og “d Po Dy ) fee 
oan ape ¥2 
34 = ” ere oe NE rn 
ie ( Ox J ay ) (G7) 
4 ; Xp! om | ¥2 
——— = eed 
1 ( ay pee ro NE (G8) 
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So, substituting terms in the expression for average power, 


Rue (44) = 5 Re| J) (Gx Hy - Gy A) da 
£ Ql ({{lie je 


(4) a oge* 
-(je = sor — Piet Se ie. eeu )paxay] 


Taking the « , ae Fourier transform of the average power 


(4) (4) * (2) % 
op... , 22 
By IES ey ox } oe 





expression, 


(4) (x) * 2 
Reve, if] ae Re| If{( («eB" (uu mi@re = ripe - +nwe® uy (G10) 
(2) 


(h) (e\* 
ms (iP E oe Kw Ve ® cue,4)\- “0 Buy “joe 2 yt dady] 


where the Fourier transform of 2 /ay has been applied, i.e., 


F [3 ]- Ser. (G11) 


From Neth (G10), performing the operations a 
Rae 4) y= [I “Bee Dey ron (a4) + aap, BY Nong) OF a 4) | 
& * | 
~ wep 4 cm pop 2S ane “@ Bhs, 4) a (G12) 
~ jock = oe (ey) + jsp PF aiey + jake OM, y) 22 aay 
As far as the first four terms in the above expression for 
average power are concerned, they will be real. However, for the 
remaining four terms, since they correspond to cross-power terms, 
a more careful analysis has to be done. A detailed exposure to such 


analysis is given in Appendix G. 


The real part of the average power turns out to be, 
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eek? Y% (| Beale hey = Be | cil =) 


+sinh 2%y [oe (e_ Btu) CMe) + pro Bey | cM ecy| x 2)|t lady : 


and similarly, for ¥, imaginary, 
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It is clear from the above expressions that one can integrate out 
the y-dependence of the average power, leaving only the e¢ -dependence. 
The integration details are shown in Appendix G, the resulting 


expressions being, 


+%& 


ue (ox\ = am) {op Cz yoh)e fe. any are 5%) 
+ pe([AM al" [>“a!?)| 
2 
+1( oe 2 lee (Bm) + Ca) ) + 0 ( [BM arf? + Jo MG 2) sinh 2p d 


ee 
+ A («*+ 1) fe, (<x) a BN) +e ( feiGalee Boy |? )} 


(G14) 
+ (each 2% d~1)[e, Ba CMa +p, att -Y, ] 
+ ee, ( (BM Car] cee) - Bea) [oMe)( \( k?-p ) a } a} 
and, similarly, for ve imaginary, 
ee“ 2 oe a {ye ¥- ve. (0a, ze Du) (Gl4a) 
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Fig. ll. Z - direction surface current density 
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+p. (JAMaal® + 1D“ Gal?) | 


ce = 
+8 (S +%"Jfer (CMe + BMH) + 0. (Gea? +1 BMC |?) | cin 2%" 
2 





Pap? f 
2 oa eae | 2\7 
Bal (x Ys le (ca) ~ Bx) + Ds (1c etf2-f Be] )3] 


rt 2 
- Zuy¥," (Bec cs] -| BM oa] CCl) ke\ d \ dx 
where in both equations (G14) and (Gl4a), the equality %G = ¥, ; 


has been applied. 


H. THE CHARACTERISTIC IMPEDANCE IN TERMS OF THE 
DISPIERSION CHARACTERISTICS 


In general, one could express the average power in terms of the 
current in the coplanar strips as, 


{ 2 
Pits Pte (H1) 


from where, 


1 Que (H2) 


ie i 


a 


In the present work, it was assumed in the dispersion character- 
istic part that the surface currents through the strips were uniform; 


furthermore, the current can be expressed, as shown in Fig. 11, as 
LebMeal =f (Ux,-2x,) dx 
strip 


=f Iai] dm & W 
strip 


where the approximation is valid since uniform surface current 


(H3) 


densities were assumed. 
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Furthermore, since it was also assumed a unit amplitude square 


pulse as the axial surface current density, i.e., J (x) = 1, it follows 


zm 
that, 
, | (H4) 
= lj | 
Te. we | (H4a) 


All constants in the average power expressions of equations 
(G14) and (Gl4a) are completely defined by choosing appropriate 
values for the width of the strips, separation between strips, thick- 
ness of the slab, dielectric's permittivity, and frequency of operation. 
The propagation constant pe is now well-defined from the value 
of the corresponding effective dielectric's wavelength » lew 


2z 
 . (H5) 


where » is determined from the dispersion characteristic curves of 
mis. /, 8, and 9, 

Therefore, the expressions to be evaluated for CPS characteristic 
impedance are, 


Re ae f 1“ (e-wfe (No) + D%e) 


+o (IAS(oa}? + 1D“ cea]?)] 


+d (Z-y Ye) ] ey (Bead + <M CM) 2p. (18 \uy|*\oag[?) |sinh 2%ed HG) 
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2 
+(eosh 2¥2d -t\}<, BEN my CO (ul + Pe JBM ey[[c“ce]] ( 2 = “| 


4 Ser, (1B ea] clay - By fo Maa (L2 ~ 2) df dee 
and, similarly, for the case where ; is imaginary, 
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+3 (pe tH") [ee (CMa + Bd) +120 [CMa] +1 BM Gul?)] sin 2"a 


Ed (n2_¥," 2) [ee (COG) _ Buy) rs je (ie amet eal) 
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I. NUMERICAL INTEGRATION 

The numerical integration problem, as far as the characteristic 
impedance is concerned, is the evaluation of equations (H6) and (H6a). 
Again, although the limits of integration in eo are ~- co and to, a 
prior knowledge of the behaviour of the integrands allows the use of 
truncated limits. 

As in the dispersion characteristic part, the numerical integra- 
tion method used is a modified Simpson's rule. Practical truncation 
limits are again -~1700 and +1700 in steps of & = 0.5, which will give 
6800 computing points for a given set of parameters. 

Once the physical parameters are specified, namely, the sub- 
strate's relative permittivity, the separation between conductors to 
width of the conductors ratio, and the substrate's thickness to free- 
Space wavelength ratio, all the constants are FOR Bmad except for 


the propagation constant (3 . However, this value can be interpolated 
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from the dispersion characteristic curves shownin part E. There- 
fore, numerically-wise, the calculation of the characteristic imped- 
ance is simpler than that for the dispersion characteristics since no 


iteration procedure is required. 


J. COMPUTER PROGRAMMING AND RESULTS 
1. Computer program organization 
The characteristic impedance program, writtenin FORT- 
RAN IV language, accepts the following data: 
= Dielectric's relative permittivity 
a Width of the conductors, in millimeters 
= Thickness of the substrate, in millimeters 


- Ratio of the thickness of the substrate to free-space 
wavelength 


= Ratio of separation between conductors to width of 
conductors 


2 Ratio of effective dielectric's wavelength to free-space 
wavelength 


Likewise, oe limits of integration, and the step of integra- 
tion eee to be specified. 

The program starts by calculating all constants and para- 
meters and directly evaluates the impedance integral; basically, the 
procedure is the same as in the dispersion characteristic part except 
for the iteration search for X since now it is completely defined. 


The program is outlined in detail in Appendix H. 
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@. Computer results 

Characteristic impedance information was obtained for three 
different relative permittivities, namely, for TS IZ = lonsand 20). 

As inthe dispersion characteristic part, this information is 
presented in terms of ratios of the different parameters, i.e., S/W 
and D/>X , in order to permit the use of any set of actual physical 
parameters. 

For each permittivity, curves were obtained for values of 
S/W equal to %, 3, 1, 3/2, 2, 3, and 4, corresponding to values of 
iy from 0.01 to 0.06. 

Characteristic impedance curves for € a 12 are presented 


in Fig. 12, for € = 16 in Fig. 13, and for ©, = 20 in Fig. 14. 
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Characteristic impedance curves for €= i 


Fig. 12. 
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Fig. . Characteristic impedance curves for €= 16 
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iV. APPELCA TIONS 


The present method is applicable to problems based on any 
coplanar strip-type configuration and also to problems related to 
slot lines, suchas coupled slots or coplanar waveguides. 

For the latter case, the configuration is, in general, as shown 
in Fig. 15, and the main difference with the coplanar strip problem 
is the knowledge of the electric fields instead of the surface current 
densities. 

Referring back to equations (B46) and (B47), these can be 


expressed as, 


E. (0x ) Mo (0 2) 

7a Bx(%) My («,B) 1) 
Mi (WB) M2 («, 8) 
Ma Ce, >) My (mx, ) 


where, 
Ms (%,@) = My (se, 2) My (, (2) ~ Ma (%, 8) Ma (eB (L2) 
therefor e, 


My (8) Ma (~,2) (L3) 
ae _ Pla : x} 
Me («, 8) Ss (x) Me (~, 2) & Ce a ( 


Similarly, 
M, (s, A) Ea (ex) 
Ms (x, 6) ex (re) 


My (%,6) M2 (mB) 
Ma (“,) My (%,8) 


Jx (e) = (4) 
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and, 


- Ma (%@) M, (v<, a) _ 
Ms (x, 6) Sete Ms(x,8) Ex (a) = Ja ; (45) 


In general, for any slot-type configuration, z-directed electric 
fields can be neglected and the transverse, or x-directed electric 
fields can be approximated with square pulse-type distributions. 

Unlike the coplanar-strip problem, here it is possible to obtain 
an odd and an even type of electric field excitations across the slots. 


Ineither case, equations (L3) and (L5) reduce to, respectively, 


a Mo (x, 6) 5 
Ms (x, Bx; (oe) = 4a) (L6) 


and 


M, («, B) 


where i =e, o depending on the type of excitation, even or odd, as 


—— 


shown in Fig. 16. 


¢ 


As far as the odd case of excitation is concerned, & (~lassumes 


eeftorm similar to J (% ) in equation (Cll), i.e., 
JU © ‘éK& 6 Xx 
8x, («) = > sin > (stw)sin SW (L8) 


Performing the necessary Fourier transform and algebraic 


reductions becomes, for the even case of excitation, 


q x : 
os (x) = SS (S+w) SINS wW, (L9) 
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Applying Galerkin's method and by Parseval's theorem, in the 
same approach as for the coplanar strip problem, the final iteration 
equation for the effective dielectric's wavelength becomes, considering 


only equation (L7), 


+ O¢ x r 
i, Bx, (x) Bx, (K1 du =o (L10) 
~ ee ee | 
or, 
. Mi ) a 
J BEB beta Grd de -c way 


where M,(% ; fe) through M,( x, @) are defined in equations (122) 
through (125), respectively. 


Remaining procedure is correspondingly similar to that for the 


coplanar strip case. 
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Even and odd transverse electric fields 


58 





V. CONCLUSIONS 


The method used in this present work succeeded in providing an 
accurate set of results for the dispersion characteristic and the 
characteristic impedance of a Coplanar Strip transmission line. 

It was also discussed that an extension of this method will allow 
corresponding results for the Coupled Slots or Coplanar Waveguide 
configuration. However, for the latter case, care must be taken in 
adapting the equations to the even case of excitation due to the fact 
that the Fourier transform of the transverse electric field will now 
be real in character, instead of imaginary as in the odd case of 
excitation. 

The amount of algebraic work can be reduced substantially by 
the use of matrices, which for the present configuration will have 
dimension eight. However, the algebraic approach used here gives 
much more insight to the problem and also reduces sharply the 
required computer time. 

As far as the theoretical development is concerned, the only 
approximation used in the representation of the axial surface current 
density, i.e., a square pulse, and the value of zero assumed for the 
transverse surface current density. Depending on the desired accu- 
racy, more complete, and indeed complex, representation of these 


surface current densities can be used. 


a) 
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In the numerical integration part of the problem, the integrating 
method itself has been proven to be sufficiently accurate. In fact, 
the only significant approximations used are the choice of the trun- 
cation limits, the integration step and the iteration step. However, 
several different sets of these choices were tried and, at least, for 
the present configuration, negligible variations were observed. 

It is clear that the two programs used, i.e., the dispersion 
characteristic and the characteristic impedance, could be lumped 
together into one program reducing in about half the amount of compu- 


tation and, therefore, the required computer time. 
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APPENDIX A 


PUES yY VECTOR POTENTIAL FUNCTIONS 


Pee GONMIiC HEREZIAN VECTOR POTENTIAL FUNCTION 
Consider a homogeneous, source-free, isotropic region; hence, 
there is no charge density and 
V-6=0. (1) 
From Vector Algebra, it is known that the divergence of a vector 


is zero if the vector is, in turn, the curl of another vector; there- 


fore, one can state, 


9-9x1 =O tC) (2) 
G=-jop Vx, | (3) 


For time-varying fields 
Vx = jwe a (4) 
so applying equation (3), 
Vxih = jwe (- pop Vx Ty ) 
RGR | (5) 
eh? Wx (Ty +k? 7h) 
mY, = ie Ty, +U¢ | 
where, 
ie wpe 


Similarly, for time-varying fields, 
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9x & = ~jopr (7) 
so, applying equation (3), 
Bx (-jep DT) «~jupp (ET +d) a 


Ue. ie, — 7? Tey = he Ty =F Cp 
(9) 


Up to this point D has not been defined; arbitrarily choose 


D= Voi (10) 
Then, equation (9) becomes, 
V2 ty + AF kh =0 sl 
Also, it is known that, 
eB =O VeH# =O 
ic 27 ry, => 
PY (12 iL, + gd | (12) 
LUO. + 9-Vb =0 
Summarizing, for TM modes, 
=~ 1a) te 
ies aes v= 10h (13) 
"‘“ae 12. So 
Hah UH, + VV- icp (14) 
= Ux on rer . 
Similarly, for TE modes, 
4, 2» 7 
H= swe Vx Ip, (15) 
ie — 35.7 
= Ou «Tl 
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B. £2 AND WPRMODES FROM VECTOR POTENTIALS 


i. TE Modes 
&= “J@p 7x 1h 
i 





U=k2 Ths VV. 


OK Te 
where ( satisfies: 


> _ 2) 
roe dee TOO . 





For TE modes, & as 0; therefore, 
h h 
rol 
Bap (S-S) 
ox oy 
from where, 
b. ab 
fy = Ihe il ; 
‘Let: . 


= ~Y¥2—-> 
ic" = dP 2 Le e 
Then, from equation (19), 
KG b"+ e pr =O 
where, 


he = ye y- 


Summarizing, 


M.=k2 i (x,4) otte 
Me etye?”? & ob 


& =i Gy x 


td 
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(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 
(26) 


(27) 








2. TM Modes 
&= 42 Te +d. Re 
= Uxux ie 
% = we x ice 
For TM Modes, # 


— ~y¥2 = 
IG = G8 (x,4) 2 Wak 


~ 


and, 
via pe ad ee p< = © 
where, 
oe Ys 
Summarizing, 
& «© Ls he ent 
Getyrhn pet? 


~~ jwem = 
w= 5 Yy Cage 





= 0; therefore, 
Z 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 
(34) 


(35) 


Equations (Al) and (AZ) correspond to equations (33) and 


(25), respectively. 
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TRANSVERSE ELECTRIC AND MAGNETIC FIELDS 


Consider Maxwell's point form of Ampere's law, 


Ux = joe & (36) 
= /d%e_ 3% > sekx eke OXy  O5Sx 
rey FF) ( om OK =) +e Co oy } 
(37) 
= jwe(& A,+ Sy Gy 4 8» Ga ) 
from where, 
 OS2 - 32 
“Bq ~ oe) de Ee _ 
Ox f ote (39) 
Se gigs li 7 
Oy Oxy (40) 
a! og "iiee S 
and, due to the z-dependence of the fields, i.e., e™ , One wean 
apply, 
= Gi” wee (41) 
Sr aueine from equations (38) through (40), 
Oe 
Y iy - = — = = jwe Gy | vt 
9) of 
= _ = jwe Ga 
(44 
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Similarly, starting from Maxwell's point form of Lenz' law, 


Ux =-jop % (45) 
is o€ 
dex (Be - Sat) +24 (SE “oe Stee (St- Se) ae; 


=~ pop (Huda + hy Gy +2 Ga) 


from where, 


47 
~¥ By =~ wp (47) 


aE 
Y&- “Se jap hy (48) 
ra) ob, 
S84. S =“jopee 
of 


i) 


Substituting the value of & from equation (48) into equation 
x 


(42), one obtains, 


ase Jue 08a 3 
By YY S Va eP*y) 

2 _y Oe. ° Ga (50) 
Hey Vsaipe) = ¥ SO -j* Se 


Similarly, substituting the values of & 7 from equation (47) 
into equation (43), of He from equation (43) into equation (47), and 
of 3 from equation (42) into equation (48), one obtains, 

us 


- For # 


’ 
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oat ~ Se 


. O82. 
m IP? % (Sy + jap %x| 


Hy (¥*+w%pe) =) one + we 


{ Ove 
Ax = oa (Y Ox ty by 


_ For & ; 
y 


o G2 


2 ae 


ve 


' Ox2 
nog ¥ 8 = jops ( SS +jwe Ex) 
oste 


By (poten be)\= 


Bye (VSG 
= For <a 


OG2 
y= &, -—~Y se 


a 


Oz 
ay wer “ok 


S,. (v4 ope) = 4 oS ie, - 


OK2 


°f 


I oO 
mene (7 2 ~jwp 7m) 


summarizing, 
2 o% 
x= oe (v= “eP a 
(ss) Se es = Sy + jw wp) 252 | 
(<6) $y zz (¥ one + | a2) 
(S7) Hy os (y a a we O82 | 


Equations (A4) through (A7) correspond to equations (54) 


threugh (57), 


Peawectiy cl. 


Sy oe, NY se 
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= —jep (39 ~jve &.) 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 





ee eNO 


DERIVATION OF FIELD EQUATIONS 


The differential equation defining the field functions is, as 


expressed in equation (B5), 


< B; (wy) ~ Ve? B (wy) =O» bet,2,3, | (58) 


f 


In turn, the propagation constant ¥. is expressed as, 


(59) 


2 z 
vt api Le a id ( he; : 
As stated before, the propagation constant for air medium, 
regions 1 and 3, is always real. However, for the dielectric medium, 
the propagation constant will be imaginary for small values of & and 
Peal for large values of &. 
For region 1, the field function differential equation solution is 
as follows, 
3? 2 
og? D, (e, 4) a " ©, (6,4) =O 
(s?-¥,7) O, (s) =f) (60) 
. +Y, ~¥ 
. D (ox, 4 | =K, e ee Kee i , 
Clearly, for infinite positive value of y, the field should vanish, 


Or, 


—& ce D, (0,4) = O 


1 


] 
a ®, (x, y) =, oy on 
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ppmoiarly, for region 3, 
x YG 
oy? ®, (cc, ) =e Ds (% 4) =O (62) 


+ Ys — ¥. 
®, (%,4) = Ks @ ee a iad 
In this case, for infinite negative values of y, the field should 


vanish, or, 


Lxe D. (, )=o 
72 | (63) 


2, By(wy) =k OF, 


Finally, for region 2, there are two solutions, corresponding to 
the real or imaginary character of the propagation constant; the solu- 
tion corresponding to the imaginary case will first be sought as 


follows, 
at 2 
3 Bz (%,4) + Yo O, (x,y) = O 


(st Y,*) 0, (S)=O 
Setiy, 


(64) 


-. D, (x,y) = Ks SiN Yey +Kg cos ry , 
r , 
For values of eo greater than ea , the propagation constant 


wee be real, or 


®, («, 4) = K> sinh ey + Ke cosh Mead (65) 
where, : 

Ky = t Ks (66) 

Koike , (67) 


One can define the following constants, 
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A(x) = Ke 
B(x)= Ks 
B, («) = Ky 
acc = Ke 
D Ck) = Ka 


(68) 


then, finally, applying the functions to the electric and magnetic field 


cases, 
Oo (wy) = AN) an (69) 
om (KY) t Be sinh YY -b Cx) cosh ey (70) 
BP" (ees) = BC) Sin aq = a" (ey CaS ‘aes (70a) 
BP (uy)= Dw e ?F (71) 
DO Coc, ye Nc ef . 
io) 


¢h) _ oh) : (4) at 
DB, (uy)= B(x) sinh ey + C'(«) cosh to (73a) 


Ch) _a@ (i) 
Dr Cug}= By (x) sin YY + C (Kw) cos %Y (74) 


Y. 
B,4 wy)=2 Dem 27 
Equations (69) through (74) correspond to equations (B10) through 


(B15), respectively. | 
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fee te NDI Xx D 


MATRIX CURRENT AND FIELD EQUATIONS DERIVATION 


From equations (B32) and (B33), one could express the constants 


poe ic ) and A! a in terms of cl? 6) and c! ) x) as, 
(el kez \2 f ple 6s¥3 (kez \? _- 
A‘ (x) = (2) {¢ (| (=) sinh Y2 d + cosh od 
al 
-~C eR alte) ~L] sinh tea } (75) 


| 2 | 
Nalen a {ce X ¥ (2) a (ee yrs | sinh Yod 
wp (Z 


oe | 


(76) 
(h) “ 
Oey | (er)? [( pe LV sez + pete) sinh % 4 
hee2 \2 
+ Vo Vs (=) cosh Yo J} 
Define the following constants: 
be os [ Keo \2 
N, = (=) | 8 (z gi sinh Yo d+cosh Ye d | 
(77) 
hee \2 « Y L 2 
> = i) sev | (¥) -1] sinh Yod 
. (78) 
_ ey hee 2 €\s hee \2 = 
alt aa Ez Vz (2-4 Ratialy (79) 
ce l | 5 kez \2 ( , 
Nu wOpeY { [r) iC ll ere a wpa¥e | sinh Ye ie 
+ w. Ys th cosh Ye d } 
Therefore, equations (75) and (76) become, 
I (x) = W), CAR) is No asl (xj so 
pli (®%) = 3 Cn) ~ Wy a (hl (rx } (82) 


(al 





from where one could solve for e : )( *« ) and oc! h ) &)interms of 


Ne )( x ) and AS ; Mx) as follows, 


of (4) (83) 
8) 2 i onmeatn AD 
eae . ‘ 
h 
a! Lao ee [ws GV ae a“ cu)] (84) 
N,Nu-NeNs e 


Once again define the following constants: 


Nu 


P - —_—__— (85) 
Ni Ng -Nz2N3 

pee (86) 
Ni Ny -N2Na 

ae Na (87) 
Ni Ny- ans 

an ey (88) 
Ni Nu -N2Na 


Therefore, equations (83) and (84) can be expressed as, 


cl (en) « e NEV) Pa alil ( (89) 


which correspond to equations (B42) and (B43), respectively. 
Substituting equations (89) and (90) into equations (B34) and (B35), 
one obtains, 


(h) ee Y. . 
Reo {he (Sy (GE) 1 8 + (BEY A ) sinh 





the? Fy cosh Ye d + ke? } (91) 
F 





‘< Me { (a EY Roe (GS) B) sinned 


+t, cosh ad \ - Sy Coe] 
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— 


—— —_ 








Sa & @i 





- L 


jwen'g Sh). Sieh -—™, 


omell | eet * 








ANd {sinh %a [vey ? ~ J a ere (ey- | 2 





. - hece 
~jx¥ = (2) PR, J tcosh Hd [jwers(c) 2 (92) 
~ Inve, — joy] ( (=. 1} | -j«rf 


+ Aa) feinh %A fj oie (es 
co 


=) 1] 84 jure ( VR 
~jver Ve Ff [+cosh Yd [ Jo R ~jwe. Ye (2 JF P 
+ \nY = )-1 2} ~ SweeY, } = Ja (x) 


Define the following constants: 


Qye ke? (5 [ey] R+ 2 (Eye) sinh ka 








Wo Wy Yo z \ lcs 
(93) 
+ kee? P y Osh Yed + hee ? 
Qz + kez ( apt (=) : P-iRae (2) B ) sinh %d (94) 
+ R cosh Yed } 
Qs = j | sioh anes fe - =e igmapere rn 
hec2 \2 
~“Y 7 ¥ (z=) P| (95) 
.+ cosh Y;d | wens (sz) 13 - KY fy ~ ec¥ (2 y-1J2} 
-«y} 
Qy = | { sinh Y2 dA (i yo: -1|? + xY 2 ce 
- WE, % ® | 
+cosh ed [xy R -wee Ye ( $2 iy Ree [( )’- J 1} 3] (96) 


- wey, } 


(fe 






ie e- — 

mii 

So! ne 

—- 

@ 1 Sle Q  ———— - _ 
i- 


—— 











Therefore, equations (91) and (92) become, 
QAM ed = Q, ACC 2 JL Co (97) 


Qs AM comer + Gy AEC = Ja (ee) (98) 
e ) 


from where one could obtain expressions for a! (*&) and As M x ) 


in terms of JA ) and i) as follows, 


1 


3) = ——_-___ 
Q2, Ry * Qe Qs 


La Ja ~ Qs Jy (oo | (99) 


l 


—>—_— 
®, Qy + Q2Q23 


[az Jae) + Gy Ja (68) | (100) 


which correspond to equations (B44) and (B45), respectively. 
Substituting equations (99) and (100) into equations (B36) and (B37), 


and solving for & a «x ) and & ( )interms of as ) and Je ), one 
x 


obtains: 
2 
———'_— | Q, Ja(u) - Qz J, (ue) | = 
a az, | Je ~ Qs Ino) } = Ee) (101) 
\ 
_ ae Ja (%) -~ Q3 J 
Q,Qy + Pz Qa {i la 2 (* « ot (102) 
+ [po Y, [Q. Jac) + Qy Jo) } } a Ey (x) 
Define the following constants: 
Mix ke? Qu 
' Q:Qyt+Q2Q3 ee) 
—_ ~ he? Qs 
Q,Qy +Q2Q3 es 
Beta ien-nra) 
Q) Qy + Q2 Qa is . tae 
a ( wie Y, Qy + «¥ | (106) 
Q) Qy+Q2Qz3 . 








Then, finally, equations (101) and (102) become, 
M, (x, 2) Ja (we) + M2 (,8) Jx Ce) = &, («) (107) 
ve (%, 2) Je (&) + Mg (%) 8) Ix (od = Sy Coe} (108) 


which correspond to equations (B46) and (B47), respectively. 


2s 





/Vlegledcil 1B 0,Ga a 


Cee aN EE Ep eFOUATIONS RE-DEFINITION 


From equation (A3), the character of the propagation constant 


was defined as lossless, i.e., 


¥= 356 


So, equations (77) through (80) are modified as, 


to (EF LSE (Eefenk ia veh 


ej) Lae: {| sinhy d 


a ik (=y ss (Z)- | sinh Yed 





| (% * hoe 
Ng * @ Pot {| opete 7 mane M dees kes }P- 1] sinh % 


+ wpe Yo ( 2) cosh re df 


Consequently, equations (85) through (88) become, 


ee 
NiNg + lel N)s\ 


IN2| ; 
- \), Na+ [Nzi[Nal 


=; ie 
NiNu + [N2]Ns| 


ie) 
Ni Ny 4 | Well Val 
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(2019) 


(110) 


iy 


ears 


(ES) 


(114) 


(155) 


(116) 


(117) 





Similaniy, “quations (93) throwgh (96) are re-defined as, 


Qk {2 (BYR ape SF (ge) -1]1BL} ssh no 


+ hoz Ty cosh Yod + he? 
aaj ke? {( SE (GE) V+ HE (Ea) BI) sian o 


+ {P| cosh Ye A i 
Qazcinh Yoel {ue 2 (ZS) Reece Srl es) )-}18l f 
~cosh Yea fests (ts 3) (Pel a Py -oee % E (GEA) f + 0/2 
By = | {sinh Xd — e ae e2\2 2] ena 2 1a des) 
+ W62%, PB Jreah Ya [wall seoee% (Ls) ® 
+068 (E)EIB] + oeev } 


Finally, equations (103) through (106) become, 





et ko? Q 

© I la@.|@s~ai[oul 
M3 as | «poy, |@a| - “8 Su) 
Mu" iia aqaay MRO ePe ¥ lea | 


(118) 


(119) 


(120) 


(121) 


(122) 


(123) 


(124) 


( 125) 


Furthermore, since from equation (Cll) J (®%) is imaginary and 
Z, 


J (%)is zero, the expressions for the different constants A through 


x 


D need to be re-defined. 


A 





Equations (B44) and (B45) become, 


RG J a (ox) 
oa? oe nal 
a... Ja| 5 (cx) 
‘ = Te Jas -eilaal " ea 


7 


“ach correspond to equations (Clo) and (CI7), respectively. 
Similarly, equations (B42) and (B43) become, 


(28) 
ON ny = PA Cx) + (P, as (o)| 


which correspond, in turn, to equations (C18) and (C19), respectively. 


Finally, equations (B38) through (B41) become, 


BY! CK) = iy, [2% [es \2 ocx) +a a{( Xe 12 PVC Meus] f (130) 

Brad =; a { wa [(E Jt] Cee) + e9pe 6 (= 22 F [ceil } (131) 
(e) am kece \2 (e) 

D «) = (Ze/ Cc) (132) 
Gacy = J ( (my [oc (133) 


which correspond to equations (C20) through (C23), respectively. 
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APPENDIX G 


CHARACTERISTIC IMPEDANCE RELATED PARAMETERS 


It is convenient to recall the Hertzian potential functions for the 


three regions of analysis. The electric field functions are, 


Bony) = AM & HI 


(134) 
h fe) “ (€) P ¢ él ec) (, y. 
2, (%,Y4)= Bo) sinh Tey + CW) cosh Ye (135) 
BS (wy) = Dm @ 77 (136) 


from where the corresponding y-derivatives can be obtained as, 











(e} ~Y¥, (y-d) 
oO YM ae |! | ee 
an * ) (e} ° 
=z =, | Ba) eosh Hy + Cu) sinh rey ] (138) 
E 
\ (2) Y. 
Seely, Ge =a (139) 
oy . 


As stated in the. dispersion characteristic part of the present 
work, specifically in equations (C16) through (C23), the magnetic 
field related constants have imaginary character. Therefore, it is 


convenient to modify the magnetic field Hertzian potential functions 


to account for this fact. 


(h) . ~¥,(4-} 
Dd x,y] =) [A (| 2 1 (140) 
os (a4) = | | (a Ga sinh Yey + [cS] cosh rey | (141) 
OM oy)= j [Deal e 4 (142) 
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and, again, the corresponding y-derivatives can be obtained as, 


Ls oe _a\ 
= ( () 
—— =} “ (B\«)} cosh YY +c“ eal sinh vey | (144) 
obs (145) 


Since an integration in the y-domain comprises three well.defined 
regions, one should separate the y-integral into three complementary 
integrals, one per each region. 

Therefore, analyzing each term in the average power expression 


of equation (Gl2), one obtains, 





B94) bs 1Y) = IAM eall pale - (146) 
Buy) BO Gay) [Beal sinh? ey +1C “ml cosh Vey 
+2 BIMa) CHa) sinh Yay cosh Yoy ae) 
3 
Bi, y) “ («,4) =[Daale>*f (148) 
a Y, S 
DN (u,4) 58 (xq) = [BS eal? sinh Vey +| Meal? cosh Yey 
) BN GallCM Gal sinh voy cosh YoY (150) 
5 (x,y) Ge, u)= | iD", Fe 2%5y 
Gis) 
_ = oo" = ¥7] Aly, i a 
(152) 
30," od, a (e), {2 2 (el, {2 2 
se =X [[Bel® cosh Yey + [c'a]* sinh* Vey 
_— (153) 


eo Ce) sinh vey cosh we | 
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3b ct ripe alte 1 
"Oy OY 
ob,“ of" i y2 
ca ey 
ob a5 
oy 34 


Pare) ge (y-< 


= Y,° { (Ber? cosh” na eC me” sinh® ny 


+2 (BO calle“ Gal sinh ea cé3ah vant ] 
od, (4) ab 4 = y2 | ae |? Bae 
oy “84 


gi" -2Y; (yd) 
10,4) 22 28 HY, Aa [AM] 4 


(8) 


1 By) Se —< =% [ BMealcMrfsinh? eq + [Ba] Cee cash ty 


+ (Ba [BM Gal +00 [CMGal) sinh Yey cash You | 


Ch) & 
' (é) Y: 
J b- (m4) = ae po aaa e* 34 


¥,(y~d) 
J 2a 5, y) = =k wy [Aa] 2° ss 


ae Be, 4) =~Ye [[Boalc“lan cosh’ Ky + Ba [cM Ga sink YY 
+ (BG) (04 cen] + Cae (CM Ger] } sinh ny cash hy | 


j ~ 203° BN, Y) = =i5, D(a) | Soe ) ef 


i(y-d 


( 2 
ge ie €) Bi (0, y) = ve [ Bea]c™ («)| cosh Vey | BO Gea] cM) eink ey 
+ (BOB cea] + Coe) [CM oy) sinh Yey cosh Wey | 


ob 
4 oy 


ch) 


. 
BN w,y) = ¥y Da) DM FO 
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(154) 


(ieSie) 


(156) 


(157) 


(158) 


Eo) 


(160) 


(G1) 


(162) 


(163) 


(164) 


ClGS} 


(166) 








(e) 
yO 4) =. = Y, Au LA Coal @ ~ 2¥;(4~d] 
(67) 


bw, wh ae ~ ie [Baa foal cast rey +1 Baal c fiat, 4 
(168) 
+ (Bud [Bea] + coc Gas] ) sinh ¥ 24 cosh YoY } 


| ae 4) i ~ V5 Duy [Dy] 2? 27 (169) 


Therefore, the expression for average power in equation (Gl2) 
becomes, 


~-2¥, (y-a) 


eg) = oe ff pend P(e lanl. LAMeeP) 
r +p (* ane \(e. IDO rg]? + pe Ip“\cu))? Je Sil 
sinh? vey [seu (<2 | BM)" +p. | BM awl?) - wa Ve" (€2 [clon] *a polo («al?) 
Hoh Ye ({ Bel] Ce (ay — BEV [CU Gecsl}) + 8c KF Yo (Be Loar] (Merah!) 1 29) 
sesh’ Ye 4 [x*ew (€2 [ci*leal*+ pe Jo Gaal?) = cops ¥Q2 (ex [BI al? + poo | BB Coa}? ) 
$e Yo (Bele) [cua -[ Bay fo lued) + oc kez Yo. (EBM Ga] Ce -B ele“) 


sinh 2Yey [wp (€, BN) CC) + Pe [BO cro Ci] ) (02 ¥) 1 } a = 


and, for the case where Y, is imaginary, the following relations hold, 
(2 = Yo 
sinh ey = | sin ‘ey 
sh YY 2 cos Y, 4 
mnh* Yeu =-sin’ Y, 
oh 4 = Cos" = 
inh zy = | sin a cal 
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Therefore, equaticn (170) becomes, 


- 
Roe (mi) = se Sf {oe (#@vAN(c, Ac]? 4p. jaMeol? a 4-4 
“4 ws («*- Yn \e, [DSIca)|?4+ Vo {th carl?) e°"3/ 
~sirt y; 4 [x2 wo (<2 [BIceal* 4 pe] Bal’) 4 were (<, 1“Gr)*apelc*bal] 
- “Be xy’ ({a"a| C010) — BO Goa Jol) (170a) 
- he x4 (BEd fo Gal -L BM Gal cGa) | 
a Ye y [oc pw (<, [c“Gal*+p. \c“oal*) +a ¥"7 (€21 Bal 4p. Emo 
— oe Y,* (BM Ga |CMeeal - | BM cea] Cen) 


= he? ¥F (JBM Gal cea - Bo Mua) | dudy, 


These two expressions correspond to equations (G13) and (Gl3a), 
mespectively. 
The y-dependence of the average power disappears through 


corresponding integrations on regions l, 2, and 3, as indicated 


below: 
- Region 1, 
” ~2Y,(4-d) | 
fe dy = 3y, (171) 
5 , 
- Region 2, 
i . { 
J sinh hey Ay ==, sinh 2%d-gd (172) 
° 2 
d 2 i: Bl 
f cosh hy dy =o sinh *hdt+ad (172) 
° 2 
d 
. { 
J sinh hy dy exp (cosh 2% 4-1) (174) 
d le 
J sin Ky ty =dd- ae sin Me (175) 
° PA 
d & a 
[east ly aye bd ach, sin 2A (176) 
A é 
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d (177) 
{sin 2% dy = 1 (1-008 2¥%3'd) 


a? 


- Region 3, 


2% _ 
Je Ay = Oy, 


~®O 


(178) 


Substituting equations (171) through (178) into equations (170) and 


(170a), one obtains equations (G14) and (Gl4a). 
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